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ABSTRACT 

We present the first results of an XMM-Newton EPIC observation of tlie luminous 
X-ray source population in the face-on supergiant spiral galaxy MlOl. We have studied 
the spectral and temporal properties of the fourteen most luminous sources, all of 
which have intrinsic X-ray luminosities exceeding the Eddington limit for a I.AMq 
neutron star, with a subset in the ultraluminous X-ray source (ULX) regime (Lx 5^ 
Xo39 erg s~^). Eleven sources show evidence of short-term variability, and most vary 
by a factor of ~2-4 over a baseline of 11-24 years, providing strong evidence that 
these sources are accreting X-ray binary (XRB) systems. Our results demonstrate that 
these sources are a heterogeneous population, showing a variety of spectral shapes. 
Interestingly, there is no apparent spectral distinction between those sources above 
and below the ULX luminosity threshold. Nine sources are well-fit with either simple 
absorbed disc blackbody or powerlaw models. However in three of the four sources 
best-fit with powerlaw models, we cannot exclude the disc blackbody fits and therefore 
conclude that, coupled with their high luminosities, eight out of nine single-component 
sources are possibly high state XRBs. The nuclear source (XMM-10) has the only 
unambiguous powerlaw spectrum (F ^ 2.3), which may be evidence for the presence 
of a low-luminosity AGN (LLAGN). The remaining five sources require at least two- 
component spectral fits, with an underlying hard component that can be modelled by 
a powerlaw continuum or, in three CclSGS, Oil hot disc blackbody (Ti„=0.9-1.5keV), plus 
a soft component modelled as a cool blackbody/disc blackbody /thermal plasma. We 
have compared the spectral shapes of nine sources covered by both this observation 
and an archival 100 ks Chandra observation of MlOl; eight show behaviour typical of 
Galactic XRBs (i.e. softening with increasing luminosity), the only exception being 
a transient source (XMM-2) which shows little change in spectral hardness despite 
a factor of ~30 increase in luminosity. We find no definitive spectral signatures to 
indicate that these sources contain neutron star primaries, and conclude that they 
are likely to be stellar-mass black hole XRBs (BHXBs), with black hole masses of 
~2-23M0 if accreting at the Eddington limit. 
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1 INTRODUCTION 

The X-ray emission of spiral galaxies arises in a combination 
of discrete source populations, hot di ffuse gas and an active 
galactic nucleus (AGN) if present fe.g. lFabbianollOSsI) . How- 
ever, for many spiral galaxies, the overall energy output in 
the X-ray band is dominated by the e mission from the few 
brightest extra-nuclear point sources jRoberts fc Warwick! 
[2000 ) . The "top-heavy" X-ray luminosity functions derived 



recently from Chandra observations of nearby spiral galax- 
ies with heightened levels of star form ation activity (e.g. 
iKilgard et all l2002l : IColbert et alll2003tl serves to empha- 
sise this point. The nature of these luminous X-ray sources, 
which may individually outshine the rest of the galaxy, is 
therefore of great interest. 

The brightest of these discrete sources are the ultralu- 
minous X-ray sources (ULXs), defined as those extra- nuclear 
X-ray sources possessing X-ray luminosities > 10^^ erg s~^, 
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which exceeds the Eddington luminosity of a IAMq neu- 
tron star by a factor > 5. The first examples of this 
class were detected by Einstein (e.g. iFabbiano fc Trinchieril 
Il987l) an d they were later detected i n large numbers by 
ROSAT JColbert fc Mushotzkvl [Tool iRoberts fc Warwick! 
[2000). ASCA studies of ULXs demonstrated that they dis- 
play "high" state multicolour disc blackbody spectra, and 
in some cases spectral transitions to a "low" powerlaw- 
dominated state, features in common with many Galac- 
tic black hole X-ray binary (XRB) systems. However, the 
masses inferred for the compact objects in these ULXs, as- 
suming Eddington-limited systems, are ~ lOOM©, much 
larger than the masses of known Galactic black hole can- 
didates (c.f. iMcGlintock fc R,emillardll200.'j) . Several models 
have now been put forward to explain the high X-ray lumi- 
nosities of ULXs. These include: sub-Eddington rate accre- 
tion onto a new class of "intermediate-mass" 
black holes (IM BHs, e.g. iGolbert fc Mushotzkvl Il999t 
iKaaret et al XRBs w ith thin super-Eddington accre- 

tion disks llBegelmanI l20n2h : XRBs emitting anisotropically 
l|King et alJl2nnih : and XRBs as scaled down versions of 
quasars (mic roquasars) with relatiyistic j ets directly in our 
line of sight (|Mirabel fc Rodn'gxiedlTflQfll ). A small fraction 
of ULXs are also kno wn to be associated with recent super - 
novae (e.g. SN 1979. Ilmmler. Pietsch fc Aschen bach"l99^. 
The observational evidence does not currently exclude any of 
these models completely. However it does now suggest that 
there are two separate populations of ULXs. Comparatively 
large numbers of ULXs have been found to be associated 
with regions of current star formation i n starburst galax- 
ies (e.g.. T he Antennae (NG C 4038/391. IZezas et alJl2002l: 
NGC 3256, iLira et al.1 120021: NGC 4485/90. iRoberts et alJ 
I2OO2I: Arp299 (NGC 3690V IZezas. Ward, fc Murravll200,tl " 
and hence must be associated with the young stellar pop- 
ulations residing in these regions. Other ULXs have been 
found in elliptical/SO gala xies where the stellar popula- 
tions are much old er (e.g. Ilrwin. Athev fc Breemanl l2003t 
IColbert et alJ2003t l , although it is not entirely clear whether 
many of these sources are associated with globular clusters 
within the galaxies, or are background AGN. 

The detailed behavioural characteristics of XRBs in our 
own Galaxy may be the key to our understanding of both 
the nature of ULXs and other very luminous X-ray sources 
in nearby galaxies, particularly since at least three Galactic 
black-ho le XRBs have been seen to re ach ULX-like lumi- 
nosities jMcChntock fc RemiUardll2003ll . However, we have 
previously been unable to derive sufficient detail from the 
study of X-ray sources in nearby galaxies for the neces- 
sary comparisons, due to the photon-limited nature of the 
available data. This has changed with the large effective 
area spectro-imaging capabilities of XMM-Newton, where 
we are now able to begin probing the detailed character- 
istics of luminous X- ray sources in nearby galaxies outside 
the Local Group fe.g.lPietsch et al.ll200ltlMiller et alJl2003t 
iPietsch. Haberl fc Vogle jl2003t) . 

One particularly good target for such a study is the 
nearby supergiant spiral galaxy MlOl (NGC 5457). Its face- 
on aspect and low line-of-sight Galactic hydrogen column 
density provide an ideal opportunity to study the discrete X- 
ray source population in a galaxy similar to the Milky Way, 
away from the obscuring hydrogen in our Galactic plane. 
MlOl has previously been observed with all of the major X- 



ray observatories. Its X-ray emissio n was first detected with 
the Einst ein X-ray Observa tory (McCai nmon fc Sanders! 
Qj^84; Trinchieri. Fabbiano fc Romai ne 199C|). -RQg^T stud- 
ies revealed the presence of numerous discrete sources 
JWang. Immler fc Piet's chlll999l) as well as a substantial dif- 
fuse component (Snowd en fc Pietschlll995h . Most recently, 
the superb subarcsecond imaging resolution of Chandra has 
revealed > 100 discrete sourc e s in the central ~ 8 arcmin- 
utes alone JPence et alJ I2OOII : iMukai et aT]l2003l) . as weU 
as diffuse emission tracing the spiral arms described by a 
two-temperature thermal plasma (Kuntz ct al. 2003). This 
galaxy was also thought to host five hypernova remnants 
(HNRs), defined as luminous discrete X-ray source s with 
supernova remnant (SNR) counterparts iWanelll999tl . How- 
ever, the superior spatial resolution and astrometry of Chan- 
dra has now shown that two of these X-ray sources are not 
in fact directly coincident with the SNRs, and that another 
shows temporal vari ability indicating an accreting XRB 
JSnowden et alJl200J) . The two remaining sources (MF37 fc 
NGC5471B) were not covered by the original Chandra ob- 
servation but are within the field of view of the XMM- 
Newton observation presented in this paper. 

This is the first of a series of papers on the XMM- 
Newton observation of MlOl. Here we present the results 
of the spectral and timing properties of the luminous X- 
ray source population, restricting ourselves to those sources 
with sufficient counts to perform meaningful spectral fitting. 
All of these sources have luminosities > 2 x 10^^ erg , 
i.e. above the Eddington limit for accretion on to a I.4M0 
neutron star. We also compliment this data set with an anal- 
ysis of two archival Chandra observations of MlOl (see sec- 
tion l2.2t . The layout of this paper is as follows. In sections |5| 
& El we outline the data reduction and source selection cri- 
teria. In sections |K|we describe the spectral fitting and 
timing analysis techniques. This is followed in section |S|by 
a detailed discussion of the spectral and temporal results for 
each source, and in section |7|we discuss discuss the over- 
all properties of this bright X-ray source population. Our 
conclusions are summarized in section |S1 



2 OBSERVATIONS & DATA REDUCTION 

2.1 The XMM-Newton Observation 

MlOl was observed with XMM-Newton for 42.8 ks on the 
4th June 2002 (ObsID 0104260101) as part of the XMM- 
Newton Survey Science Centre (SSC) Consortium Guaran- 
teed Time (GT) programme. The EPIC MOS-1, MOS-2 fc 
PN cameras were operated with medium filters in "Prime 
Full Window" mode, which utilizes the full ~ 30 arcminute 
field of view of XMM-Newton. This completely encompasses 
the g alaxy's D25 isophotal ellipse diameter of 23.8 armin- 
utes llTullvlll988D and corresponds to a spati al coverage of 
62.8 k pc at the distance of MlOl (7.2 Mpc, IStetson et all 
[1998?) . The data were pipeline-processed using the SAS {Sci- 
ence Analysis Software) v5.3.2 and all subsequent data anal- 
ysis was carried out using SAS v5.4.1. Full-field light curves 
were accumulated for the three exposures to check for high 
background intervals of soft proton fiares. There are numer- 
ous small flares throughout the exposure, and the four most 
prominent peaks (each with count rates greater than ~100 
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Figure 1. XMM 3-colour image of the MlOl field (stacked MOS & PN). Colours correspond to 0.3-0.5 keV (red), 0.5-2.0 keV (green) 
& 2.0-4.5 keV (blue). The 14 luminous X-ray sources are marked with circles denoting the spectral extraction radius used. Two bright 
foreground stars in the field are marked with squares. 



counts per second in the PN) were cut from subsequent data 
analysis, leaving a net good time for each camera of 36.7 ks 
(MOS-1 & MOS-2) and 36.6 ks (PN). All data (images and 
spectra) were created using patterns corresponding to single 
& double pixel events for the PN (0 & 1-4), and patterns 
0-12 for the MOS cameras. We also set FLAG=0 to reject 
events from bad pixels and events too close to the CCD chip 
edges. 

2.2 The Chandra Observations 

In order to study the changing spectral and temporal prop- 
erties of these sources, we compliment our XMM analysis 
with recent Chandra data of MlOl. Two Chandra obser- 



vations of MlOl have been performed at the time of this 
analysis; a long (~100ks) observation on the 26th March 
2000 (ObsID 934) and a shorter (~10ks) observation on the 
29th October 2000 (ObsID 2065). T hough an analysis o f the 
longer observation is presented in iPence et al.l ll200ll) . we 
have reanalyzed this data to utilize more up to date calibra- 
tion files available in the Chandra CALDB v2.11, including 
new gain maps, improved plate scale, and improved PHA 
bin calculation. The standard data processing was done at 
the Chandra X-ray Center, and subsequent analyses were 
performed using the CIAO software {Chandra Interactive 
Analysis of Observations) v2.2.1. The data were screened 
for periods of high background leaving 98.2 ks (ObsID 934) 
and 9.6 ks (ObsID 2065) of good time. A more detailed de- 
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script ion of the data anal ysis is available in iKilgard et alJ 
l|2003h . In addition, whilst ' Pence et al.l J200ll) only include 
the back-illuminated ACIS-S3 chip in their analysis (which 
covers ~11 percent of the D25 ellipse of MlOl), we include 
a set of 5 ACIS chips for both observations; in addition to 
the S3 chip, we use the front-illuminated ACIS 12, 13, S2, 
and S4 chips. The long observation covers 51 percent of the 
D25 ellipse of MlOl, but the combined data set gives 78 per- 
cent coverage and enables us to perform spectral and timing 
analyses of discrete sources in a uniform, consistent fashion. 



3 SOURCE IDENTIFICATIONS 

In this paper, we are studying the brightest point sources 
in MlOl at the time of the XMM observation. To define our 
subset, we set a lower limit of 300 counts accumulated in the 
PN camera (0.3-10 keV), which we judge as the minimum 
requirement to perform detailed spectral and temporal anal- 
ysis. To search for such sources, we ran the SAS script EDE- 
TECTCHAIN on 0.3-10 keV cleaned PN and MOS images 
with a 3(7 detection significance threshold. We detected four- 
teen sources with >300 counts in the PN data, and we name 
them 'XMM-n', numbered as n in descending PN count rate 
order. The sources are listed with their coordinates and PN 
count rates in Tabled The only exception to the above is 
XMM-14, a bright source which was not identified in the 
PN data by the source detection algorithm due to its close 
proximity to a chip gap. It was detected in the MOS data, 
so we have derived a count rate consistent with those of the 
other sources (i.e. corrected for chip gap losses) by faking a 
PN spectrum in XSPEC using the best fit absorbed power- 
law model from its combined spectra MOS and PN spectra 
(see section llj, the response matrix and ancillary response 
files and exposure time from a nearby source (XMM-9) . The 
complete source list for this observation will be presented in 
Paper II [m preparation). 

Figure shows an XMM 3-colour image of the MlOl 
data from the three EPIC cameras. Each source is labelled 
with its 'XMM-n' source number, and the circles represent 
the spectral extraction radii used. The field of view of the 
combined Chandra dataset covers the positions of all of the 
XMM sources except the brightest source in the XMM obser- 
vation (XMM-1). A Chandra source list was created from the 
merged data using the wavdet ect algorithm in CIAO , and the 
full source list is presented in lKilgard et alJ (jioO^ together 
with comprehensive details of the source detection tech- 
niques used. Eleven of the XMM sources covered by these 
data were unresolved, while the remaining two (XMM-2 & 
10) are resolved into two discrete sources each. The official 
Chandra names (following the Chandra naming convention) 
are listed in Tabled The astrometry of Chandra is generally 
good to ~1 arcsecond on the aimpoint ACIS S-3 chip, but 
the accuracy does decrease with increasing off-axis angle. In 
general, the XMM and Chandra positions agree to within 2 
arcseconds. We have also cross-identified the positions of the 
fourteen sources with ot her X-ray observat ions of MlOl, i.e. 
ROSAT ma and PSPC JWang et alJll999l) and the original 
Chandra ACIS-S3 analysis bv lPence et alJ (1200111 . for which 
the relevant ID numbers are also listed in TableQ Note that 
XMM- 2 is coincident wit h the ROSAT PSPC source P17, 
which dWang et alJll999ll associate with HRI sources H22 



& H25. However, the XMM position closely corresponds to 
the position of H25 (offset ~7 arcseconds) rather than H22 
(offset ^^33 arcseconds). 

Figure |5| shows the source positions overlaid on an op- 
tical blue DSS2 image to illustrate their locations in MlOl. 
iMatonick k. FesenI ([1993) identified 93 SNRs in MlOl, aU 
of which are within the XMM field of view. We have listed 
two of these (denoted MF) which we currently believe to 
be co-located wit h our X-ray sources. However, the Chan- 
dra observations dSnowden et aljboOlf) have demonstrated 
that two other SNRs, MF54 & MF57, originally associated 
with the ROSAT Wl sources H29 (XMM-9) & H30 (XMM- 
14), are in fact too far offset from the SNR positions to be 
Eissociated with them. Three of our sources are associated 
with the five giant HII complexes in MlOl (first detected 
in X-rays with ROSAT by Williams_& Chu 1995), strongly 
suggesting that they are in some way associated with the 
star formation activity occurring within them. 



4 SPECTRAL PROPERTIES 

The XMM X-ray spectra for each source were extracted 
using circular regions enclosing all of the visible emission, 
of radii varying between 10 and 35 arcseconds. The small- 
est apertures were used in order to avoid contamination 
from nearby sources in crowded regions of the galaxy, or for 
sources located next to a chip gap (e.g. XMM-9, 10 & 14). 
The background regions for each source were taken using as 
large an area as possible, close to each source without con- 
tamination from other sources and with approximately the 
same DET-Y distance as the source region (to ensure simi- 
lar low-energy noise subtraction). We used the SAS task ES- 
PECGET to simultaneously extract source and background 
spectra and to create response matrices (RMFs) and ancil- 
lary response files (ARFs) for each source. The spectra were 
then binned to a minimum of 20 counts per bin in order to 
optimise the data for statistics. 

The spectral analysis of the XMM data in this study has 
been performed using XSPEC vll.1.0. All errors are given 
at the 90 percent confidence level unless stated otherwise. 
The MOS and PN spectra have been fitted simultaneously 
in the 0.3-10 keV band for each source. We included a free 
normalisation constant to account for differences in the flux 
calibration of the three EPIC cameras, finding the relative 
normalisations to be typically within ~ 5 percent. 

We initially fitted the spectra with simple absorbed sin- 
gle component spectral models: powerlaw continuum (PL), 
MEKAL optically-thin thermal plasma, conventional black 
body (BB) and a multicolour disc black body (DISKBB). 
This last mo del describes an accr e tion disc in its high /sof t 
(HS) state jMitsuda et alj Il984l : iMakishima et alJ Tl98^ . 
Each model had both a free absorption component plus an 
additional one fixed at the line of sight Galactic hydroge n 
column density (1.16x10^° cm~^. lDickev fc Lockmanll990f) . 
In Table l^l we show the best fit parameters for the PL con- 
tinuum and DISKBB models for each source, together with 
the observed fluxes^ and observed plus unabsorbed lumi- 
nosities (0.3-10 keV) from the best fit model (highlighted in 

^ The fluxes quoted here are an average of those measured in the 
three EPIC cameras, except for instances where the normalisa- 
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Figure 2. XMM source positions (circles) overlaid on a DSS2 optical blue image (note that source circles do not represent positional 
error radii and smaller circles are used to separate close sources). Large square and circular footprints denote the fields of view of the 
PN and MOS cameras respectively. 



Table 1. Source positions and cross- identifications. 



Source 


RA 


Doc 




PN Count Rate 


Rosat ID" 




Chandra ID 


Other ID^''' 




(J2000) 


(J2000) 


(cts ks-i) 




Pence'' 




This work 




XMM-1 


14''03™03.9= 


+54° 


27' 


33" 


182.90 


P13/H19 








MF37 


XMM-2 


14''03"n4.3* 


+54° 


18' 


05" 


138.00 


P17/H25 


P48 


CXOU 


J140313.9+541811 


















P51 


CXOU 


J140314.3+541807 




XMM-3 


U*" 04'" 14. 2* 


+54° 


26' 


03" 


128.84 


P28/H45 




CXOU 


J140414.3+542604 




XMM-4 


U*" 02'" 28. 5* 


+54° 


16' 


24" 


75.47 


P8/H10 




CXOU 


J140228.3+541626 


NGC 5447 


XMM-5 


14''04'"16.7* 


+54° 


16' 


13" 


36.66 


P30/H47 




CXOU 


J140416.8+541615 




XMM-6 


14'' 02'" 22. 4^ 


+54° 


17' 


58" 


32.85 


P6/H9 




CXOU 


J140222.2+541756 




XMM-7 


U** 04'" 29.1^ 


+54° 


23' 


53" 


28.01 


P32/H49 




CXOU 


J140429.1+542353 


NGC 5471 


XMM-8 


14''02"'03.5* 


+54° 


18' 


28" 


27.49 


P4/H3 




CXOU 


J140203.6+541830 




XMM-9 


U*" 03^21. 5* 


+54° 


19' 


46" 


23.95 


P19/H29 


P70 


CXOU 


J140321. 5+541946 




XMM-10 


U*" 03^12. 4^ 


+54° 


20' 


55" 


21.06 


P16/H23 


P38 


CXOU 


J140312.5+542053 


















P40 


CXOU 


J140312.5+542057 




XMM-11 


U*" 03"' 36.0* 


+54° 


19' 


24" 


19.86 


P21/H36 


P104 


CXOU 


J140336.0+541925 


MF83 


XMM-12 


U*" 02™ 52. 8* 


+54° 


21' 


10" 


17.60 


P11/H18 


P5 


CXOU 


J140252.9+542112 




XMM-13 


M*" 03^41. 2* 


+54° 


19' 


03" 


12.11 


P22/H37 


P107 


CXOU 


J140341. 3+541904 


NGC 5461 


XMM-14 


U*" 03™ 24.0* 


+54° 


19' 


49" 


10.57 


P19/H30 


P76 


CXOU 


J140324.2+541949 





Notes: °- ROSAT ID numbers correspond to PSPC (P) & HRI (H) detections CWang et al. 19991 . ''Chandra ID numbers from 
Pence et al. (2001). '^'MF' indicates SNR detection of Matonick fc Fesen (1997). "^'NGC' denotes giant HIT regions e.g. Williams & 
Chu (1995). 
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Table 2. XMM-Newton single-component model fits. 



Source PI? DISKBB" fP 









K 




r 






T„, (kcV) 




Obs 


Obs 


Unabs 


XMM- 


■1 


2, 


c--|+l).24 


1, 


■^"-0.07 


305.4/291 




-I no + U.UV 

-^••^•^-0.07 


290.1/291 


5.01 


31.1 


34.2 


XMM- 


■2 


0, 


ycr + O.lS 
' 15 


2, 


1 0+0.09 
— 08 


262.5/235* 


[~0.12]^- 


[~0.64]'= 


600.9/235 


3.64 


22.6 


28.5 


XMM- 


■3 


1, 


7t;+0:32 
■ ' '■'—0 30 


2, 


vo+0:i9 
■ '°-0 17 


168.1/155* 


<0.17 


0.58+q'q3 


185.4/155 


2.34 


14.5 


31.5 


XMM- 


■4 


2, 


Q4+o:4i 

-0.39 


2, 


2o+o:i4 

■^'-'-0.06 


156.5/151 


<0.33 


J^-UJ^-0.07 


143.6/151 


2.13 


13.2 


13.5 


XMM- 


■5 


1, 


47+0.79 
—0 82 


2, 


Qg+0.38 
— 27 


37.0/25* 


<0.59 


0-98^02? 


48.0/25 


1.11 


6.9 


9.7 


XMM- 


■6 


1, 


oR+0.57 
^0-0.59 


2, 


t:c:+0.28 
■OO_0.24 


35.0/38 


<0.31 


79+0.09 
"• '^-0.08 


40.4/38 


0.68 


4.2 


8.1 


XMM- 


■7 


5, 


79+4.01 
'^-2.26 


7, 


97+2.73 
■^'-1.76 


66.8/54 




'^•^'-'-0.04 


52.1/54 


0.40 


2.5 


10.1 


XMM- 


■8 


2, 


1 c: + 0.86 
^"^-0.93 


2. 


nfi+0.34 
■^"-0.29 


44.1/33 


<0.79 


1 oi+°-2° 

^•"^-0.21 


43.9/33* 


0.61 


3.8 


3.9 


XMM- 


■9 


3, 


77+0.94 
' '-0.80 


2, 


■^"-0.21 


34.3/31 


1 99+0.63 
-^-■^^-0.52 


96+°-i5 

'-'•^"-0.13 


30.6/31 


0.57 


3.5 


4.2 


XMM- 


■10 




<0.81 


2, 


9t;+0.41 
■^''-0.17 


17.9/19 


[~0.12]^- 


[~0.45f 


39.2/19 


0.53 


3.3 


3.4 


XMM- 


■11 


1, 


fic+o.gg 

"0-0.51 


3, 


Qc+0.55 
"^"-0.35 


64.0/44* 


<0.72 


34+0-05 
'-'••^^-0.06 


69.0/44 


0.33 


2.1 


6.2 


XMM- 


■12 


2, 


"^-1.02 


1, 


on+0.23 
■^'J-0.25 


44.1/48 


90+° ®° 


1 45+0-49 
-^•^"^-0.28 


48.8/48 


0.53 


3.3 


4.5 


XMM- 


■13 


8, 


40+306 
^'■'-2.65 




>6.16 


64.6/37 


o 17+2.17 
■^■^'-1.01 


16+0-03 
'-'•^'-'-0.04 


56.3/37* 


0.21 


1.3 


10.1 


XMM- 


■14 


4, 


,0 + 1.75 
J^^-1.21 


1, 


of.+0.27 
■S'3-0.23 


10.0/22 


1 fio + 1.04 


1 40+0.35 
^•^^-0.26 


14.5/22 


1.06 


6.6 


9.1 



Notes: "Spectral models: PL=powerlaw continuum model and DISKBB=multicolour disc black body model. ''Observed 
flux in the 0.3-10 keV band, in units of 10^^"^ erg cm~2 s~^. '^Observed and unabsorbed luminosities in the 0.3— lOkeV 
band, in units of 10^* erg (assuming a distance to MlOl of 7.2 Mpc). ''Total absorption column (including 
Galactic), in units of 10^^ cm~^. "^Parameter unconstrained due to high xi- *x2 > 1 therefore fits with two- 
component models are attempted (see Table. The best-fit model for each source is highlighted in bold. 



Table 3. XMM-Newton two-component model fits 

Source ModeP fcT/Ti„ f X^Jdoi Ax^" l-P(F-test)= Flux Ratic? 



[keV) Soft Hard 



XMM-2 


PL-I-DISKBB 


67+°-''^ 
"•"'-0.19 


0. 


qn+0-U4 

'■'"-0.06 


1. 


41+0.25 
'^-'--0.18 


199. 


,4/233 


63.1 


100.0 percent 


0. 


,26 


0. 


,74 




PL-I-BBODY 


0.38l°:^f, 


0. 


9-1 +0.03 
^^-0.03 


1, 


CI +0.07 
■'^-'--0.18 


200 


.7/233 


61.8 


100.0 percent 


0, 


.20 


0, 


,80 




PL-I-MEKAL 


47+0.15 
-0.15 


0. 


70+0.11 
'°-0.18 


1, 


Qq+0.09 
■='•'-0.08 


240 


.5/233 


22.0 


>99.9 percent 


0, 


.05 


0, 


,95 


XMM-3 


PL4-MEKAL 


1 o~+0.34 
■'■•"''-0.34 


0. 


gg + 0.33 

'"'-0.27 


2. 


5Q+0.21 


160. 


,3/153 


7.8 


97.4 percent 


0. 


,05 


0. 


,95 




PL+BBODY 


1 97+0.51 
^■^'-0.46 


0. 


cQ+0.23 
"''-0.21 


3, 


1 0+0. 36 
■-'-•'-0.29 


164 


.6/153 


3.5 


79.5 percent 


0, 


,17 


0, 


,83 




PL-I-DISKBB 


1.941?,:?! 


0. 


oQ+0.47 
''y-0.34 


3, 


99+1.10 
■^^-0.60 


165 


.0/153 


3.1 


76.1 percent 


0, 


.27 


0, 


,73 


XMM-5 


PL-I-MEKAL 


6.54t^-» 


0. 


1 7+0.15 
'-^'-0.05 


2. 


jO+0.47 
■-^-"-0.44 


26. 


,5/23 


10.5 


97.8 percent 


0. 


,20 


0. 


,80 




PL-I-BBODY 


1 01+4.79 
-^-■^-^-0.64 


0. 


17+0.12 
-^'-0.08 


1, 


40+0.87 
■^•'-1.05 


30 


.5/23 


6.5 


89.1 percent 


0, 


.15 


0, 


,85 




PL-fDISKBB 


n 90 + 3.49 
■'-^='-2.62 


0, 


1 7+0.30 
' -0.06 


1, 


R1 +0.62 
■"-^-1.40 


30 


.8/23 


6.2 


87.8 percent 


0, 


.16 


0, 


,84 


XMM-11 


PL-I-BBODY 


<1.19 


0. 


1Q+0.03 
'-^'^-0.04 


1. 


7.0 + 1.01 
■ '"'-l.OO 


55. 


,6/42 


8.4 


94.8 percent 


0. 


,44 


0. 


,56 




PL+DISKBB 


0.83+°?° 


0. 


90+0.09 
■^■'-0.06 


1, 


-,0 + 1.08 
■68_j 24 


56 


.5/42 


7.5 


92.9 percent 


0, 


,49 


0, 


,51 




PL+MEKAL 


o.78ij;j° 


0. 


99+0. 32 
■^='-0.25 


2, 


gi+0.87 

■^-'--0.32 


56 


.8/42 


7.2 


91.9 percent 


0, 


,15 


0, 


,85 


XMM-13 


PL-I-BBODY 


3.4lt2-35 


0. 


1 9+0.03 
'-'-^-0.03 


0. 


of5+1.19 
■■"'-1.56 


34. 


,3/35 


30.3 


>99.9 percent 


0. 


,40 


0. 


,60 




PL-fDISKBB 




0. 


1 0+0. 04 
■^''-0.03 


0, 


g^ + 1.02 
■"^-1.88 


35 


.4/35 


29.2 


>99.9 percent 


0, 


,41 


0, 


,59 




PL-I-MEKAL 


p. 70+0.76 

^•''^-1.81 


0. 


1 0+0. 06 
■-^"-0.03 


0, 


ocr + 1.22 
■''0-1.33 


38 


.4/35 


26.2 


>99.9 percent 


0, 


,43 


0, 


,57 



Notes: "Spectral models and parameters as in Table 01 except : BBODY=blackbody model and MEKAL=thermal plasma (solar 
abundances). ''Improvement in the statistic over the single component model powerlaw (PL) fit, for two extra degrees of 
freedom. '^F-test statistical probability of improvement of fit over single component model (PL). Fraction of total flux in the 
soft and hard (PL) model components over 0.3— lOkeV band. The best-fit model for each source is highlighted in bold. 



bold). We omit the BB fits from the table as they did not 
provide significantly better fits than the DfSKBB model, 
except in the cases of XMM-7 and XMM-13, -wliich we -will 
discuss individually in sections 16.71 and 16.131 We also omit 
the MEKAL fits from the table because, in most cases, tlie 
spectra -were too hard to be fitted well by this model alone; 
two exceptions to this are XMM-7 and XMM-14, both of 
which will be discussed individually in section |5] 

Tiie spectra of eight of the fourteen sources were ade- 
quately fit with tiiese simple absorbed one-component mod- 
els with X? i$ li fi-nd tlieir spectra plus best-fit model are 

tion constant of the PN differed from the MOS cameras by >10 
percent, in which case the average MOS flux is quoted. 



plotted in Figures |3&|3] together with the ratios from botfi 
PL and DISKBB models, fn the remaining cases wiiere the 
xi > 1, we went on to fit the spectra with more complex 
two-component models. The results are shown in Tablej^for 
the five sources where a significantly better fit at ^ 95 per- 
cent confidence level was achieved^, with the best-fit model 
highlighted in bold. Figure|^shows the spectra/ratios for the 
best-fit two-component models on the left, while the ratios 
of the data/model for the remaining fits are shown on the 
right for comparison. Where the two-component fits did not 



^ No significant (^ 95 percent) improvement was made to the fit 
of XMM-8 with the addition of a soft component. 
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Table 4. Chandra single component model fits. 





I, n a n Hvn 




PL° 










^ X 




c 

X 




CXOU ID 





r 


2 /J f 
X /dot 


Nh 


Tin (keV) 


2 /J f 
X /dof 


Obs 


Obs 


Unabs 


XMM-1 






















XMM-2 


J140313. 9+541811 






















J 14Uol4.oH-04ioU 1 


|~1.2iJ 


r.^ ,0 /I Ol e 




|~U.i2J 


['^ U.D4J 


3/ r 


U.Uo 


U.4o 


U. I y 


j\.iVliVl-0 
VViViiVi-^ 


71/10/11 /I 'i-L^AOfifiA 
J 14U414.o-t-04ZDU4 


o c:i+0.35 
•^■^^-0.33 




114/1 "jn 


1^ C7+0.23 


'J-^'^-0.04 


1041:/ IZU 


z. / 


1 7 Ozl 


1 1 Q 71 
1 ly . / 1 


XMM-5 


J140416. 8+541615 


U.OO_g gg 


— 0.23 


30 /58 


<0.30 


^■^'-'-0.16 


38/58 


1.04 


6.43 


7.88 


XMM-6 


J140222. 2+541756 




















XMM-7 


J140429.1+542353 




















XMM-8 


J140203.6+541830 


- 


- 


- 


- 


- 


- 


- 


- 


- 


XMM-9 


J140321. 5+541946 






29/45 


54+0-38 
'-'•'^*-0.34 


1 07+018 


25/45 


0.50 


3.09 


3.40 


XMM-10 


J140312.5+542053 


<2.25 


1 7O+0.63 
'•^-0.40 


4/13 


<0.61 


1 17+0.33 
^-^'-0.36 


6/13 


0.18 


1.09 


1.29 




J140312.5+542057 


<3.16 


^•"'-0.47 


13/13 


<1.27 


53+014 


19/13 


0.13 


0.79 


0.93 


XMM-11 


J140336.0+541925 




^•°°-0.37 


24/34 


<0.24 


Q 44+0.11 
"-^^-0.05 


42/34 


0.30 


1.83 


2.92 


XMM-12 


J140252.9+542112 


'J-y'i_o.47 


^•^''-0.23 


30/36 


h0.15]<= 


hi. 61]'= 


32/36 


0.52 


3.24 


3.72 


XMM-13 


J140341. 3+541904 


, -, 37+3.34 


10 n+2-27 


33/24 


5 34+2.09 


13+0-04 
"-^■^-0.02 


17/24 


0.12 


0.72 


29.23 


XMM-14 


J140324.2+541949 


4 40+111 
^■^l'-0.95 


, 79+0.23 
l-'^-0.21 


32/53 


r, 71+0.62 
^- ' -^-0.62 


, c:c: + 0.28 
l-a3_0.08 


26/53 


0.87 


5.41 


6.67 



Notes: "Spectral models and parameters as in TableEl ^Observed flux in the 0.3-8 keV band, in units of 10~13 erg cm s . 
'^Observed and unabsorbed luminosities in the 0.3-8 keV band, in units of 10^8 erg s~i (assuming a distance to MlOl of 7.2 Mpc). 
"^Total absorption column in units of 10^1 cm"^. "^Parameter unconstrained. The best-fit model is highlighted in bold. 



Table 5. XMM-Newton short-term variability test results. 



Source Bin Size Standard deviation, a statistic K-S statistic 

(s) Expected" Observed" Po-(var) x^/dof P^2 (var) Px-s(var) 



XMM-1 


200 


2.88 


X 


10" 


2 


(3.35 ± 


0.18) 


X 


10" 


2 


99.2 percent 


206/171 


96.5 percent 




XMM-2 


200 


2.72 


X 


10" 


2 


(4.07 ± 


0.22) 


X 


10" 


2 


>99.9 percent 


427/171 


>99.9 percent 




XMM-3 


400 


1.38 


X 


10^ 


2 


(1.66 ± 


0.13) 


X 


10" 


2 


97.6 percent 


104/86 






XMM-4 


300 


1.65 


X 


10- 


2 


(2.19 ± 


0.14) 


X 


10- 


2 


>99.9 percent 


149/115 


98.2 percent 




XMM-5 


2000 


2.62 


X 


10- 


3 


(3.76 ± 


0.61) 


X 


10" 


3 


97.7 percent 


29/18 


95.2 percent 




XMM-6 


1200 


4.33 


X 


10- 


3 


(5.14 ± 


0.66) 


X 


10" 


3 




35/29 






XMM-7 


1500 


3.29 


X 


10- 


3 


(3.90 ± 


0.56) 


X 


10- 


3 




12/23 






XMM-8 


2000 


2.68 


X 


10- 


3 


(4.06 ± 


0.66) 


X 


10- 


3 


96.4 percent 


23/18 




97.4 percent 


XMM-9 


1500 


3.50 


X 


10" 


3 


(4.22 ± 


0.61) 


X 


10" 


3 




30/23 






XMM-10 


2000 


2.43 


X 


10" 


3 


(3.62 ± 


0.59) 


X 


10" 


3 


95.8 percent 


34/18 


98.7 percent 




XMM-11 


1200 


4.28 


X 


10- 


3 


(7.57 ± 


0.98) 


X 


10" 


3 


99.9 percent 


73/29 


>99.9 percent 


>99.9 percent 


XMM-12 


1500 


3.35 


X 


10- 


3 


(5.57 ± 


0.80) 


X 


10" 


3 


99.4 percent 


42/23 


>99.9 percent 




XMM-13 


1500 


3.42 


X 


10" 


3 


(4.88 ± 


0.70) 


X 


10" 


3 


96.1 percent 


26/23 




99.0 percent 


XMM-14 


2000 


2.60 


X 


10" 


3 


(3.31 ± 


0.54) 


X 


10" 


3 




22/18 




97.5 percent 



Notes: "Expected and observed standard deviations from mean count rate (count s i). We only show P(var)>95 percent to highlight 
variability. 



bring the xl value down to 1, more complex models were at- 
tempted. In only one case did this significantly improve the 
fit, although the model is partially unconstrained (XMM-13, 
see section 16.131 . The results of the spectral fitting for each 
source are discussed in section |5] 

For direct comparison of the spectral shapes of the 
sources between the Chandra and XMM observations, we 
have fitted the same set of basic spectral models to the data 
from the first (long) Chandra observation, which provided 
sufficient counts for spectral fitting. We have fit the data of 
ten Chandra sources, nine of which have > 300 counts. The 
data for CXOU J140314.3-H541807 (XMM-2) is below this 
threshold (178 counts), but we have also fit these data for 
comparison purposes (although the fit parameters are un- 
constrained). Standard CIAO routines were used to extract 
source and background spectra from the reprocessed event 
lists and to construct RMF and ARF files. The ARF files 



were also corrected for the ACTS low energy degradation us- 
ing the corrarf tool. The spectral fitting was performed in 
Sherpa using chi-squared statistics with the Gehrels error 
approximation, which returns a slightly lower value of the 
fit statistic as compared with traditional chi-squared statis- 
tics. Errors were calculated for 90 percent confidence using 
the projection function and Levenberg-Marquardt optimiza- 
tion. The fit results are shown in Table |1] and all sources 
are well fit with the single-component models. However, for 
comparison purposes, we have also fit the Chandra spectra 
of the five sources best described by more complex models 
in the XMM analysis with the same two-component mod- 
els. The only statistically significant improvement over a 
single-component fit is found in the PL-I-MEKAL fit to 
CXOU J140414.3-f542604 (XMM-3), which we discuss in 
section 16.31 We compare and discuss the XMM and Chan- 
dra spectral fits in section im 
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XMM-6 Powerlaw (r=2.55) 



XMM-10 Powerlaw (r=2.25) 





Energy (keV) 
XMM-12 Powerlaw (r=1.80) 



Energy (keV) 
XMM-14 Powerlaw (r=l.l 





Energy (keV) 



Energy (keV) 



Figure 3. XMM spectra of the powerlaw sources. PN data points and the best fit model are shown in black; the MOS data are shown 
in grey. Ratios are shown for both the PL and DISKBB fits. 



5 TIMING ANALYSIS 
5.1 Short-term variability 

The X-ray variability of each source in the XMM data was 
tested by deriving short-term light curves. The data from the 
three EPIC cameras were co-added to improve the signal- 
to-noise ratio, and the exposure times were tailored so that 
each bin had at least 20 counts after background subtraction, 
giving temporal resolutions ranging from 200-2000 seconds. 
The four time intervals of high background flaring were ex- 
cluded, and we corrected the exposure times and count rates 
in the incomplete bins adjacent to the gaps with a simple 
scaling factor, excluding bins with < 0.3 times the expo- 
sure remaining of the original bin size. The resulting light 
curves are shown in Fig. |S| [left] , with error bars correspond- 
ing to la deviations assuming Gaussian statistics. Two tests 
were applied to these light curves. Firstly, we calculated a 
standard deviation of the count rate per bin from the mean 
count rate and compared this to the expected deviation of 
±18-20 percent expected from Gaussian counting noise. Sec- 
ondly, we performed a test to search for large amplitude 



variability against the hypothesis of a constant count rate. 
Table |K| shows the results of these tests together with the 
probability that the data are variable, P(var), in the cases 
where it exceeds 95 percent. These tests on the binned data 
revealed that ten sources show some evidence of variability, 
of which four (XMM-2, 4, 11 & 12) show variability which 
is significant at the Scr confidence level (P(var)>99.73). 

In order to search for any indications of gradual small 
amplitude variations in the source count rates, we conducted 
Kolmogorov-Smirnov (K-S) tests using PN light curves with 
a time resolution of 1 second'^. For each source, we compared 
the observed background-subtracted cumulative photon ar- 
rival distribution with the expected distribution if the fiux 
was constant. Of the fourteen sources, four show evidence for 
variability, but only one (XMM-11) is statistically variable 
at the 3(T confidence level. Initially, source XMM-7 also ex- 
hibited variability in the K-S test, but further investigation 



^ Only the PN data were used for the K-S test as the MOS cam- 
eras are limited to 2.6 second time resolution in Full Window 
mode. 



© 2003 RAS, MNRAS OOO.inOTI 



XMM-1 Disc Blackbody (T|^=1.33 keV) 
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XMM-4 Disc Blackbody (T, =1.01 keV) 
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XMM-8 Disc Blackbody (Tj =1.01 keV) 
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Figure 4. XMM spectra and ratios of tlie disc blaclcbody sources. PN data points and the best fit model are sliown in black; the MOS 
data are shown in grey. 
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XMM-2 Powerlaw+Disc Blackbody (r=1.41, Tjjj=0.30 keV) 
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Figure 5. [Left] XMM spectra and ratios of best-fit two- component models. PN data points and the best fit model are shown in black; 
the MOS data are shown in grey. [Right] Ratios of all other two-component models for the corresponding source on the left. 
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XMM-13 Powerlaw+Blackbody {r=0.36, ^^=0.13 keV) 




Energy (keV) 
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Figure 5 - continued 



showed that this could be attributed to variability in the un- 
derlying background flux in this region of the detector, which 
constituted 57 percent of the total (source+background) flux 
for this source. 

Similar tests on the Chandra data detect significant 
short-term variability only in CXOU J140336.0-I-541925 
(XMM-11), with P(var)>99.9 in both the and K-S tests. 

5.2 Long-term variability 

We have taken advantage of the fact that MlOl is a well- 
studied galaxy that has been observed in the X-ray regime 
with the Einstem, ROSAT, ASCA and Chandra observato- 
ries. We have used these data (except ASCA due to its poor 
angular resolution) to construct long-term light curves for 
the fourteen sources. For consistency between the datasets, 
we have plotted the fluxes in terms of the observed 0.5-2 keV 
flux, as this is the energy range covered by all the missions. 
The XMM and Chandra fluxes are measured directly from 
the best-flt absorbed single component models. We have 
used webPIMMS to derive fluxes for the Einstem IPC and 
ROSAT PSPC/HRI observations using the observed count 
rates, assuming the absorbed PL continuum slope measured 
in the XMM observation and normalized to the 0.5-2 keV 
band . In the two Etnstem IPC observations iXrinchieri et alJ 
Il99(f) . the central region of the galaxy is not resolved, but 
emission is detected from three of the XMM sources away 
from the centre of the galaxy (XMM-1, 3 & 4). AU four- 
teen sources were detected in the ROSAT PSPC observa- 
ti on, and our fluxes are derived from the count rates given 
in lWang et"al] ( 199g^■ However, si nce the ROSAT VLRl count 
rates quoted in Wang et al.l 1^9^ are a combination of four 
individual observations, we have re-analysed the archival 
data (as per iRoberts fc Warwic k 2000) and measured count 
rates (or upper limits) for each individual data set. Finally, 
we correct the count rate of the ROSAT PSPC source P19, 
which is resolved into two sources by both the ROSAT HRI 
(H29 & H30) and XMM (XMM-9 & 14). Since the fluxes of 
these two sources are almost identical in the flrst HRI ob- 
servation, we have attributed half of the PSPC flux to each. 
Note also that as XMM-2 and XMM- 10 are resolved into 
two sources in the Chandra observations, the Chandra data 
points represent their combined flux. 

The resulting light curves are shown in Figure |S| [right]. 



and although poorly sampled with only 8-10 observations in 
11-24 years, the majority of the sources do show evidence of 
long-term variability /transient behaviour. These results are 
discussed in section 17m 



6 INDIVIDUAL SOURCE PROPERTIES 

In this section, we discuss the spectral and timing proper- 
ties of each source in turn. For consistency, we quote the 
luminosities for both the XMM and Chandra data in the 
0.3-8 keV band. An overview of the spectral and temporal 
properties for all sources is given in section 17.11 



6.1 XMM-1 

This source is spatially coincident with the SNR MF37 in an 
outer spiral arm of the galaxy. However, the ROSAT data 
showed evidence of strong variability, implying that it is a 
single accreting object, perhaps an XRB formed in the su- 
pernova explosion. It is the brightest source in the XMM ob- 
servation, and its spectral shape (disc black-body model 
with Ti„=1.33keV) plus indication of short-term variability 
strongly support the hypothesis that this is a single XRB. 
Its long-term lightcurve also shows signiflcant variability of 
a factor of ~3-4 between observations, with an increase in 
the XMM observation putting it well into the ULX regime 
(intrinsic Lx = 3.3 x 10'^^ erg s~^). Similar high inner accre- 
tion disc temperatures have been observed in other ULXs, 
e.g. the ASCA study of lMakishima et alJ J2000h . This is the 
only source that was not covered by either of the Chan- 
dra observations. 



6.2 XMM-2 

This source is located in a spiral arm region contain- 
ing numerous HII regions. It is coincident with the 
ROSAT source H25(P17), but is resolved into two sources 
(CXOU J140313.9+541811 & CXOU J140314.3-f 541807) 
in the Chandra observations. However, only CXOU 
J140314. 3-1-541807 has sufficient counts for spectral fit- 
ting and it is marginally best fit with an (unconstrained) 
absorbed PL model (r=2.4) with an unabsorbed X- 



ray luminosity of Lx 



7.9 X lO"*' erg s"\ although 
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Figure 6. [Left] XMM siiort-term light curves, obtained by summing tlie counts from the 3 EPIC cameras. The mean count rate is shown 
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detection was made. 
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a disc blackbody model gives a statistically indistin- 
guishable fit (Ti„=0.6keV). The second source (CXOU 
J140313.9+541811) has a luminosity of Lx = 3.4 x 
10^^ erg s~^ if we assume the spectral fit parameters of the 
XMM PL model described below. 

Although we expect that the emission we detect in the 
XMM observation is the unresolved flux from both Chan- 
dra sources, XMM-2 is ~30 times more luminous than ei- 
ther of the ROSAT or two Chandra sources, and the peak of 
the emission corresponds more closely to the position of the 
brighest Chandra source (CXOU J140314.3+541807, offset 
~1.8 arcseconds). The J't'MM spectrum of this source is best 
described by a two-component model (see Tablesl^l&OJ- The 
spectral fit is improved over that of a simple PL fit with the 
addition of either a DISKBB, BBODY or MEKAL compo- 
nent. The PL-I-DISKBB is statistically the best fit, but it is 
almost indistinguishable from the PL-(-BBODY fit with sim- 
ilar temperatures and PL slopes. The short-term light curve 
of this observation shows definite variability, indicating that 
it is dominated by a single accreting source. Its unabsorbed 
luminosity in this observation is Lx = 2.7 x 10'^'"* erg s~^, 
classing it as a ULX. The sharp increase in luminosity in the 
XMM data compared with the previous observations leads 
us to conclude that this is an X-ray transient source which is 
likely to be related to the active star formation in its vicinity. 



6.3 XMM-3 

Unlike the previous sources, this is not associated with a 
region of obvious star formation, but it is located in a re- 
gion in the direction of the end of a spiral arm (Fig. I^J. 
The ROSAT data showed that it was strongly variable by 
a factor of ~2 on timescales of days. The Chandra data are 
well fit with a soft PL continuum (r=3.6), which yields a 
high unabsorbed luminosity {Lx = 1.7 x 10^^ erg s~^ ob- 
served; Lx ~ 1.2 X 10*" erg s^^ instrinsic), although there 
is a large uncertainty associated with the intrinsic luminos- 
ity due to the large absorption correction {Nh = 3.5 x 10^^ 
cm~^). It is also in the ULX regime in the XMM observa- 
tion {Lx = 3.1 X 10^® erg s~^), and the spectrum is best 
described with an absorbed PL-I-MEKAL thermal plasma 
model (r=2.6, fcT=0.98keV). Even though the x^ = l-05 
for this spectral model, we tested whether the fit would 
be improved by adding an extra MEKAL component, but 
this made the fit marginally worse. However, if we fit the 
Chandra spectrum with a PL-(-MK model we get a sig- 
nificantly improved fit {NH=3.23to ll, kT^O.SOtoM^f^'^, 
r=3.42lQ jg, xi ~ 0-9) which gives additional evidence of 
the presence of thermal plasma in the vicinity of this source. 
We detect no definite variability in its short-term light curve, 
but its long-term light curve shows that it has varied by a 
factor of ~4 between observations spanning ~24 years. This 
variable behaviour shows that this source is almost certainly 
dominated by an XRB. Although the MEKAL component 
only contributes 5 percent of the total flux of this source, it 
has a high luminosity of Lx ~ 1.6 x 10"^* erg s~^. The source 
is too far away from the centre of the galaxy to be contam- 
inated by diffuse emission from the disc. Neither is it in a 
star-forming region of MlOl and no SNR has been detected 
at this position. We therefore speculate that the thermal 
plasma could originate in a hot photoionized plasma (e.g. 



stellar wind) surrounding the XRB, if the companion is a 
high- mass star (see section 17^ . 

6.4 XMM-4 

This source is located in the giant HII com plex NGC 5447 in 
a spiral arm of MlOl. IWang et alj 1^9^ originally classed 
it an XRB, based on the hard spectral shape implied by 
the ROSAT PSPC harness ratios. This source was detected 
in the short second Chandra observation, although it has 
too few counts for spectral fltting. This source has bright- 
ened considerably (factor of ~3) in the XMM observation 
{Lx = 1.3 X 10"^" erg s~^), and its spectra are best fit 
with absorbed DISKBB model with a high disc tempera- 
tu re (Tin = 1.0 keV). agai n similar to the ULXs measured 
bv lMakishima et alJ 1 I2OOOI) . Both the short-term XMM and 
long-term light curves show evidence of variability; this cou- 
pled with its spectral shape leads us to conclude that the 
emission is due to an XRB embedded in the HII region. 



6.5 XMM-5 

This source is located away from the main body of th e galaxy 
and was classed as an AGN bv lWang etall (ll999D accord- 
ing to its hardness ratios, the fact that it has a blue stellar 
counterpart and lack of short-term variability. It is however 
in the direction of an outer spiral arm and may be a source 
associated with MlOl. In the Chandra observation, its spec- 
trum is well fit with an absorbed PL (r=1.8). The best fit 
single component model to the XMM data is an absorbed PL 
with an unabsorbed luminosity of Lx = 9.1 x 10^* erg s~^, 
but the introduction a cool MEKAL component does im- 
prove the fit (r=2.1, fcT=0.17keV). The MEKAL compo- 
nent contributes a substantial amount (20 percent) to the 
total flux from this source, which could indicate the pres- 
ence of substantial star formation activity in this region, 
or might again be attributable to the source itself. Its long- 
term light curve shows a factor of ^^3 variability between the 
different observations, and there is some evidence of short- 
term variability in the XMM data. This raises the possibly 
that it is an accreting XRB in MlOl, although we cannot 
discount the possibility that this is a background AGN, as 
rapid X-ray variability on the times cales of hours is know n 
to exist some in BL Lac objects (e.g. lTaeliaferri et al .1200311 . 
According to the USNO-A2 catalogue, the apparent B mag- 
nitude of the optical counterpart is 18.7, which, assuming a 
F,^ ~ spectrum typical for an AGN, equates to an ap- 
parent R magnitude of ~17.8. Using the 2-10 keV flux of 

6.6 X 10"^* erg cm"'^ s~^, this in turn equates to an X-ray- 
to-optical flux ratio log{Fx / Fr)=-0.56, which is within the 
range measured for optically identifled AGN i n the XMM- 
Newton hard X-ray survey of IPiconcelli et all (l2003l) . How- 
ever, such objects do not typically show MEKAL compo- 
nents in their X-ray spectra. Alternatively, if the optical 
source is located in the outer spiral arm of MlOl, the im- 
plied absolute magnitude of Mb = —10.6 is consistent with 
the range measured in the V -band for stellar cluste rs in the 
Antennae (-9 > My > -14. IWhitmore et al.lll999^ . There- 
fore, although this source may be a background AGN, its 
short- and long-term variability plus the presence of ther- 
mal plasma in its vicinity makes it much more likely that 
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this is an XRB associated with star formation activity in 
MlOl. 



6.6 XMM-6 

This source is located in the same spiral arm as XMM-4 
close to an HII region. Although no short-term variability is 
detected in the ROSAT HRI data, this source showed hard 
X-ray spectral characteristics in the PSPC data which lead 
the authors to conclude that it is an XRB. It was detected 
in the short Chandra observation, with an unabsorbed lu- 
minosity of Lx = 3.3 X 10^* erg s"\ The XMM spectrum 
is best fit with a soft absorbed PL model (r=2.6) with an 
intrinsic luminosity of Lx = 7.9 x 10^* erg s~^. Its short- 
term light curve doesn't show any indication of variability, 
although its long-term light curve does. The observed fiux 
peaks in the first ROSAT HRI observation, gradually fades 
up to the Chandra observation and sharply increases in the 
JfMM observation by a factor of ~3. This variable behaviour 
leads us to conclude that this source is an XRB. 



6.7 XMM-7 

This source is one considered by I Wan 3 (^9^ to be a HNR 
candidate. The ROSAT HRI source (H49) is positionally 
coincident within la error radius of the optical position of 
NGC 5471B, one of the SNRs detected in N GC 5 4 71 by 
[gkillman,. ( 198^ in the radi o and Chu & Kcn nicutd (|l98d[ l 
in the optical. [wan3 ((l993) found that the ROSAT PSPC 
data was consistent with a Raymond-Smith thermal plasma 
with a temperature of 0.29 keV and an X-ray luminosity of 
Lx = 3 X 10^* erg s~^ (0.5-2 keV). Using these characteris- 
tics, they employed the Sedov solution to estimate the total 
thermal energy of the SNR to be ~ 10^^ ergs (if that is 
where the X-ray emission originates). This is ten times that 
expected from a typical supernova and so have suggested 
that it could be a HNR. No temporal variability was de- 
tected in the ROSAT data, consistent with the hypernova 
scenario. 

More recently, IChen et alJ l)2002l) have studied 
NGC 5471B using Hubble Space Telescope (HST) WFPC2 
contiunuum band and emission lines image together with 
echelle spectra for kinematic studies. They found three 
[SII]-enhanced shells in NGC 5471, the brightest of which 
is NGC 5471B with an expansion velocity of ^ 330 km s~^, 
a kinetic energy of 5 x 10^^ ergs and hence an explosion 
energy > lO'^^ ergs. The WFPC2 images showed that it has 
a complex environment of many concentrations of stars and 
that the shell itself encompasses two young OB associations. 
The ROSAT HRI contours are centered on NGC 5471B on 
the HST images and the authors conclude that the SNR 
within it was produced by a hypernova since a shell moving 
at that velocity is unmatched by any comparable shells 
in other galaxies that have been produced by multiple 
supernovae. 

This source is an off-axis detection in the short Chan- 
dra observation, and hence has insufficient counts for spec- 
tral fitting or an accurate determination of its position. The 
XMM spectrum shows that this source is very soft. It is 
weU fit with both an absorbed cool BBODY (ri„=0.14 keV, 
xS=0.93) model with an unabsorbed luminosity of Lx = 



5.2 X 10^* erg s~^ and a DISKBB (ri„=0.16 keV, x^=0.96) 
model with an unabsorbed luminosity of Lx = 1.0 x 
10"^^ erg s~^. This cool temperature classifies this source 
as a luminous supersoft s ource (SSS), similar to those 
found in other galaxies fe.g. lDi Stefano fc Konell2003ar) . Al- 
though not shown in Table |21 this source was also ade- 
quately fit with an absorbed single MEKAL plasma model 
(fcr =0.1 keV, x^=l-2). Its light curve shows no variability 
over the timescale of the XMM observation, although there 
is evidence of variability in the long-term light curve (factor 
of ~2.5). The nature of this source therefore remains open 
to debate. A reasonable MEKAL X-ray spectrum and a lack 
of variability are still consistent with an interpretation as a 
HNR. However, the similarity of its spectral form to other 
SSSs, suspected to be XRBs, offers an alternative solution. 



6.8 XMM-8 

XMM-8 is located away from the main body of the galaxy, 
but it is in the direction of one of the outer spiral arms 
(see Fig. HJ. It was detected in the ROSAT observations, 
but showed no variability. It was also detected in the short 
Chandra observation, although with insufficient counts for 
spectral fitting. However, the XMM spectra are best fit with 
an absorbed DISKBB {Tin = 1.01 keV) with an unabsorbed 
luminosity of Lx = 3.9 x 10^** erg s~^. Even though this fit 
has a value of 1.33, the addition of a second component 
does not significantly improve the fit. There is some evi- 
dence of variability in the short-term XMM light curve, and 
the long-term data show a factor of ~2 variability between 
observations, with a sharp increase in fiux in the XMM ob- 
servation. We therefore conclude that this source is likely to 
be an XRB. 



6.9 XMM-9 

This source was orginally identified by IWan3 l|l999^ as be- 
ing possibly coincident with SNR MF54. Its ROSAT X- 
ray position was offset from the optical position of MF54 
by 6.2 arcseconds {2a error radius) and therefore consid- 
ered to be a possible HNR due to its high X-ray luminosity 
{Lx = 1.4 X 10^* erg s~^) and lack of temporal variability. 
However, this source was detected in the lo ng Chandra ob- 
servat ion, and it was clearly demonstrated bv lSnowden et al] 
i200lf) that its position (accurate to ~0.5 arcseconds) is in- 
consistent with the optical position of MF54 with an offset 
of 8.3 arcseconds between them. They did detect a faint X- 
ray source (P67) at the position of MF54 with an X-ray 
luminosity of Lx = 4.3 x 10"^® ergs s~^ (thermal model 
with kT — 0.73) which is reasonable for a SNR with a 
diameter of ~20 pc as determined by jMatonick & Fose^ 
il997l) . Our analysis of the Chandra data for the brighter 
source show that its spectrum is best fit with a hot DISKBB 
model (ri„=l.lkeV) with an unabsorbed X-ray luminosity 
of Lx = 3.4 X 10^* erg s~^, although it is also well fit with 
an absorbed PL (r=2.0). 

The XMM spectra are consistent with this scenario, be- 
ing best fit with an DISKBB model with a similar tem- 
perature (Tin =0.96 keV) with an unabsorbed luminosity of 
Lx = 4.2 X 10^* erg s"\ Its short-term XMM light curve 
shows no indication of variability, and likewise its long-term 
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lightcurve appears fairly constant over the 11 year period. 
However, its spectral shape leads us to conclude that it is 
an accreting XRB, with an inner disc temperature similar 
to those ULXs observed with ASCA. 



6.10 XMM-10 

This is the source located at the optical nucleus of 
MlOl. In the ROSAT data it has an extended morphol- 
ogy and shows no variability. In the Chandra observa- 
tion, the nucleus is resolved into two point sources (CXOU 
J140312. 5+542053 & CXOU J140312. 5+542057), both of 
which are b e st fit with absorbed PL models (r=l. 7/2.1). 
iPence et al.l i200ll) have shown that the most northerly 
source (CXOU J140312. 5+542057) coincides with the nu- 
cleus, the other (CXOU J140312. 5+542053) coincides with 
a cluster of bright stars 3.1" (~110 pc) to the south. Neither 
source shows temporal variability and they have similar un- 
absorbed luminosties of Lx = 1.3 x 10^* erg s"^ (CXOU 
J140312.5+542053) & Lx = 9.3 x 10^^ erg s"^ (CXOU 
J140312.5+542057). 

The two sources are not resolved in the XMM data, 
but the properties of XMM-10 are fairly consistent with 
their combined emission. The spectra are best fit with an 
absorbed PL model (r=2.3) with a luminosity of Lx = 
3.3 X 10^* ergs s"\ However, the XMM short-term light 
curve does show some evidence for variability over the du- 
ration of the observation. The long-term lightcurve shows 
that the combined flux from these two sources has remained 
relatively constant over the 1996-2002 period, but there is 
some variability apparent in the early data (factor of ~2), 
indicating that at least one of the two sources is a variable 
XRB system. 

6.11 XMM-11 

This i s another of the five X-ray sources considered bv lWan j 
il999i) to be a HNR candidate based on its positional co- 
incidence with SNR MF83 and its high X-ray luminosity 
(Lx = 1.2 X 10^* ergs s~^). With the ROSAT HRI and 
PSPC data. lWand 11199911 was unable to constrain the shape 
of its X-ray spectrum. However, he comments that if the X- 
ray emission is difl^use (i.e. a thermal plasma), then its X-ray 
luminosity suggests a total thermal energy of 6x 10^^ ergs. 
The HNR scenario was supported by the lack of temporal 
variability in their data. 

MF83 is located between two spiral arms (see 
Fig. 13, and is one of the largest SNRs identified by 
[Matonic k fc FesenI in MlOl by its high [SII]/Hq ratio 

and bright [OIII] emission, clearly suggesting that the opti- 
cal nebula is shock heated. Ho wever, based on ground-based 
and HST optical observations. iLai et alJ (|200Jl showed that 
MF83 is a star-forming region consisting of a '-^270 pc ion- 
ized gas shell with an expansion velocity of ~50 km s^^. 
It has composite sources at its centre which are likely to be 
OB associations and four HII regions along its rim. The large 
size of the shell and the fact that it is centered on groups of 
stars suggest that this shell is a superbubble similar to those 
found in the Large Magellanic Cloud, rather than a SNR. 

The Chandra data have been able to constrain the X-ray 
spectrum of this source. It is best fitted by an absorbed PL 



(r=2.7), and is poorly fitted by thermal emission model. 
Significant short-term variability was found in the Chan- 
dra light curve (xS ~4.3), which rules out the hyp ernova sce- 
nario. The variability, together with the PL spectral shape 
suggests that the origin of the X-ray emission is likely to be 
an XRB associated with the superbubble. The position of 
the Chandra source is coincident with the MF source posi- 
tion to ~ 3 arcseconds. 

The XMM spectra of this source are not adequately 
fit with a simple one-component model, but the only two- 
component model which shows a significant improvement is 
the PL+BB (r=1.7, kT = 0.19 keV). The PL+DISKBB fit 
(r ~1.7, Tin = 0.22 keV) gives an improvement at the 93 
percent significance lev el, with parameters similar to IMBH 
candidates studied by iMiller et alJ Its short-term 

XMM light curve also shows significant variability, support- 
ing the idea that this is not a HNR. The long-term lightcurve 
also shows some variability (factor of ~2). It has an unab- 
sorbed luminosity of Lx = 6.2 x 10^^ erg s~^ and we con- 
clude that this source is an XRB associated with the super- 
bubble. 



6.12 XMM-12 

This source wass classed as an 'interarm' AGN by 
I Wang et alJ il999l) based on its hardness ratios, blue op- 
tical counterpart and lack of temporal variability in the 
ROSAT HRI data. The Chandra data are well fit with ei- 
ther an absorbed PL model (r=1.5) or a DISKBB model 
(Lin=1.61 keV) an d no variability is evident. However, 
iMukai et al.l i2003ll do note that the Chandra position is 
5 arcseconds fro m the optica l posit ion of the blue counter- 
part reported bv lWang et alJ lll999l) and so the source may 
indeed be located within a spiral arm in MlOl. 

The XMM data are also best fit with an absorbed 
PL model (r=1.8) with an unabsorbed luminosity of 
Lx = 4.1 x 10^® erg s"\ although the DISKBB also gives 
a very good fit to the data (Ti„ = 1.45 keV). Its long-term 
light curve appears fairly constant, but its short-term 
lightcurve shows strong evidence for variability, leading us 
to conclude that it is an XRB. 



6.13 XMM-13 

This source is coincident with NGC 5461, another of 
the bright giant HII regions in MlOl. Thermal and non- 
thermal radio emission has been d etected in NGC 5461 
l)Graeve. IGg inj_^^J)^[i^l£^inskiM]M and the kinematic 
studv b\ nOu^^Cenntoira ~ |l98m the presence 

of high- velocity (455 km s"^ FWZI) gas near its cen- 
tre. However, no SNR was detected at that position by 
iMatonick fc FesenI 1^23) ■ The presence of high- velocity gas 
could be explained by a wind-blown bubble around stars em- 
bedded in the HII region, although the velocities measured 
here are far great er than those detected in other galaxies 
such as the LMC dChu fc Kennicutt|[l98^ . 

The ROSAT observations characterized this source as 
diffuse, with no temporal variability. In the Chandra obser- 
vation, the spectrum is best fit with an absorbed cool disc 
blackbody model (ri„=0.13keV), probably indicating the 
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presence of a massive (or IMBH) black hole (although no 
short-term variability was detected) . The XMM data show 
similar results, but the spectra are not adequately fitted with 
any single component model (the best-fit single component 
model is a BB with A;T=0.14keV, xl = 1-49). The best fit 
two-component model is an absorbed PL-f-BB model with an 
extremely hard PL slope and cool BB temperature (r=0.4, 
fcT=0.f 2 keV), but a PL-f-DISKBB model gives almost iden- 
tical fit statistics with similar parameters. A fit with a 
PL-I-MEKAL component also gives an improved fit over 
the one-component model, and a PL-|-(2xMEKAL model 
(r=0.3, fcr=0. 3/1.0 keV) gives a very good fit {xl = 0.96), 
but these data could not constrain the temperature of the 
second MEKAL component. There is some indication of 
short-term variability in the XMM observation, and the 
long-term lightcurve shows that it is also varying over longer 
timescales (factor of ~5), demonstrating that this source 
is likely to be dominated by an XRB. The unabsorbed lu- 
minosity in this observation is Lx = 1.0 x 10^^ erg s~^, 
putting this source into the ULX regime. We also speculate 
that although partially unconstrained, the two temperature 
MEKAL thermal plasma components could well be emission 
from the superbubble surrounding the XRB. 



6.14 XMM-14 

This is another of the HNR candidates of IWand Jl999l) . 
based on its close proximity to the optical position of SNR 
MF57 (4.0 arcseconds offset). Its X-ray properties in the 
ROSAT observation were similar to those of XMM-9 (H29) 
with a high X-ray luminosity and lack of variability. But 
again, the accurate Chandra position of this source showed 
a clear offset of 4.8 arcseconds from the SN R position and 
is the refore discounted as a HNR candidate iSnowden et alJ 
I2OO J) . The Chandra data for this source are best fit with a 
hot disc blackbody model (Ti„ — 1.55 keV) with an unab- 
sorbed X-ray luminosity of Lx = 6.7 x 10'^* erg s~^, again 
in the range of temperatures measured in the ASCA study. 

The XMM data for this source are best fit with ab- 
sorbed PL (r = 1.86) with an unabsorbed luminosity of 
Lx ~ 8.3 X 10"^* ergs s~^, although a disc blackbody model 
also provides a good fit with a temperature consistent with 
the Chandra data {Ti„ — 1.49 keV). Although not shown in 
Table ^ it is also well fit by an absorbed MEKAL model 
(fcT = 6.4 keV, xi = 0.49), though this temperature is un- 
realistically high. The K-S test indicates possible short-term 
variability, and the long-term data show that it is definitely 
varying over long timescales (factor of ~3). We therefore 
conclude that this is another XRB. 



7 DISCUSSION 

7.1 Overview of Source Properties 

We have performed spectral and temporal analyses on the 
fourteen most luminous sources in MlOl at the time of the 
XMM-Newton observation. The most striking result of the 
spectral analysis is that the sources appear to be a hetero- 
geneous population. Five sources are best described by a 
simple absorbed DISKBB model with temperatures ranging 
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Figure 7. Plot of the ratio between the highest and low- 
est observed luminosities (Lmax/Lmin) in the 0.5— 2keV band. 
The symbols denote PL (squares), DISKBB (circles), super- 
soft (cross), PL+DISKBB/BBODY (triangles) and PL+MEKAL 
(stars) sources. 

from 0.16-1.33 keV, three of which have unabsorbed lumi- 
nosities in the ULX regime. Another four are well described 
by absorbed PL models with photon indices ranging from 
1.80-2.55. The PL spectral fits of the remaining five sources 
are improved by the addition of a second soft thermal com- 
ponent. On the other hand, eleven of the sources show some 
short-term temporal variability, and the majority show long- 
term variability /transient behaviour. These properties are 
consistent with most (if not all) of the sources in this sam- 
ple being XRBs. 

We investigate whether any trends in behaviour with 
spectral type are present in the following manner. We quan- 
tify the long-term variability by plotting the ratio between 
the maximum and minimum observed luminosities in the 
long-term data, indicating the different spectral types with 
different symbols (Figure |7|l. In each case, we have used the 
observed minimum and maximum values except for XMM- 
14, where the upper limit of the third ROSAT IfRI data 
point is the lowest value and hence the emission must be 
less than this. While two sources (XMM-9 & 12) only show 
minor fiux variability (factor of ~l-2), it is clear that the 
majority of sources vary between a factor of ~2-4, with no 
distinct grouping of sources with certain spectral types. It 
is however interesting to note that the two sources with the 
highest degree of long-term variability (the transient XMM- 
2 and XMM- 13) are those with PL-^DISKBB/BBODY com- 
posite spectra. 

We can make a simple comparison in spectral shapes 
between the single component model fits for the long Chan- 
dra observation and the XMM-Newton observation of MlOl 
(Tables |21 & EJ . The nine sources with sufficient counts for 
spectral fitting in both observations are all best fit with 
the same (PL or DISKBB) continuum model in both ob- 
servations, except XMM-14 where both models are statis- 
tically acceptable due to the low numbers of counts. Also, 
the majority of these sources show no change in powerlaw 
slope or disc temperature within the 90 percent confidence 
errors. The only significant change in shape in the whole 
sample is in XMM-3; the powerlaw slope hardens between 
Chandra and XMM observations (3.60 to 2.78) with only 
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Intrinsic X-ray Luminosity (10'^^ erg s~', 0.3-8 keV) 

Figure 8. Plot of the ratio between the hard (^2— Sfcev) soft (i^o.5— 2fcev) intrinsic (unabsorbed) luminosities vs the broad-band 
(0.3-8 keV) intrinsic luminosity for both the XMM [left] (open symbols) and the long Chandra [centre] (filled symbols) observations. 
[Right] The spectral changes are shown for those sources with spectral fits in both observations. The symbols denote PL (squares), 
DISKBB (circles), supersoft (cross), PL+DISKBB/BBODY (triangles) and PL+MEKAL (stars) sources. 



a marginal decrease in flux, but this may be offset by the 
difference in best-fit column densities of the two models. 
To illustrate the spectral hardness of these sources, in Fig- 
ure|H]we have plotted the ratio between the hard (L2-8fcev) 
and soft (I/o.3-2fcev) intrinsic (unabsorbed) luminosities ver- 
sus the broad-band (0.3-8 keV) intrinsic luminosity for both 
XMM [left] and the long Chandra [centre] observations mea- 
sured directly from the spectral fits. Figure m [right] il- 
lustrates the changes in hardness and luminosity between 
sources with spectral fits in both observations. The majority 
of sources show a softening with increasing luminosity typi- 
cal of Galactic black hole XRBs (BHXBs) such as Cyg X-1 
(e.g. McClintock & Rcmillard 2 003) . Si milar behaviour has 
also been observed in some ULXs (e.g. iKubota et alj|200lt 
iLa Parola et al.ll200ll) . though othe r ULXs are known to 
display the opposit e behaviour (e.g. iFabbiano et aljl2003al : 
iRoberts et allbooal) . The exception to this trend in MlOl is 
the transient source XMM-2, which shows little change in 
hardness despite a substantial change in overall luminosity. 

One interesting result from inspecting Figure |H1 and 
in particular the left panel, is that there does not appear 
to be any change in either the underlying spectral shape 
or in the spectral hardness across the Lx = 10^^ erg s~^ 
ULX threshold. Albeit on the basis of a small sample, this 
suggests that, in the absence of a luminosity measurement, 
ULXs are spectrally indistinguishable from somewhat less- 
luminous sources in nearby galaxies. This in turn suggests 
that sources both above and immediately below the ULX 
threshold are intrinsically very similar, if not the same, 
and that no new source class is appearing in MlOl above 
Lx = 10^^ erg s~\ Indeed, if we make the assumption 
that our sources are spherically accreting BHXBs at the 
Eddington limit (Lg = 1.5 x 1O^*(M/M0) erg s~\ e.g. 
iMakishima et aPbOOOt) . their intrinsic X-ray luminosities in 
the XMM observation imply black hole masses in the range 
of ~2-23M0*. This is reasonably consistent with the range 



of mea sured masses for Galactic stellar -mass black holes (3- 
ISAfrT^. lMcClintock fc Remillardlf2003l) . implying that there 
is no strong requirement for IMBHs to explain any of the 
present observations. It is therefore very likely that the lu- 
minous point sources we see in MlOl are simply the extreme 
end of the "ordinary" XRB population in the galaxy. This is 
consistent with the cumulative luminosity functions of com- 
pact sources in other galaxies with high star formation rates, 
which extend t o ULX luminosities and show no break at 
~ 10^^ erg s-i jGrimm. Gilfanov fc Sunvaevll2003D . 

The locations of these sources can provide clues to 
the nature of the X-ray emission. However, they can also 
be misleading if any of the sources are background AGN 
shining through MlOl. We estimate the probability of 
any of our detections being background sources using the 
hard band (2-10 keV) integral \og{N)-\og{S) relationship of 
iGampana et alJ ll200j) (we judge this to be a better indica- 
tor than the soft band log(Ai')-log(S) due to the high, spa- 
tially variable neutral absorption column associated with the 
disc of MlOl). After converting our source detection thresh- 
old of 300 counts in the 0.3-10 keV band into a 2-10 keV 
flux for a typical AGN, and folding this through both the 
exposure map for the observation (to correct for sensitiv- 
ity changes across the EPIG field-of-view) and the 2-10 keV 
log(A'^)-log(S'), we predict a total contamination of ~2 back- 
ground AGN out of fourteen sources. Hence the likely con- 
tamination of the luminous X-ray source population with 
background AGN is small. 

So are any of the fourteen sources good candidates to 
be a background AGN? The majority of the sources in 
this sample are actually firmly associated with star form- 
ing HII regions or are located within spiral arms, arguing 
that they are intrinsic to MlOl. Interestingly, the outlying 
sources that are arguably the best candidate AGN (XMM-3, 
5 & 8) are also all located in the direction of faint outer arms 
of the galaxy. The relatively high neutral hydrogen columns 



* Of course, the inferred luminosity of XMM-3 in the Chandra ob- 
servation implies a much higher black hole mass, though a large 
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correction used (see Section . 



18 L. Jenkins et al. 



we measure in the spectral fits are consistent (within the 
90 percent confidence limits) with most of the sources being 
embedded within the neutral absorbing material of MlOl, 
which ranges in column between 1 x 10^" cm~^ in the 
interarm disc regions up to ~ IQ^^ cm~^ in the spiral arms 
l|BrauJl99!j ^.^ One final test of whether a source is a good 
candidate for a background AGN is if it has a bright, point- 
like optical counterpart. This is the case for the source 
XMM-5, though its X-ray characteristics may argue against 
an AGN interpretation (c.f. Section f6. St . It is therefore dif- 
ficult to identify any of the fourteen sources as an obvious 
background AGN. 

7.2 The nature of the compact object 

If all of our fourteen sources are accreting XRB systems, 
what is the nature of the compact object? Characteristics of 
Galactic BHXBs such as state transitions and rapid tem- 
poral variability have also been found in systems known 
to contain a neutron star primary. Even though all the 
sources have intrinsic luminosities above the Eddington 
limit for accretion on to a IAMq neutron star {L%^ = 
2 X 10^* erg s~^), we know that both neutron star XRB 
(NSXB) and BHXB systems can achieve super-Eddington 
luminosities. For example, the three Galactic BHXB sources 
(V4641 Sgr, 4U 1543-47 & GRS 1915+105) have been ob- 
served with peak luminosities well into the ULX regime 
iMcClintock fc Remillardl2003^ . Super-Eddington luminosi- 
ties have also been observed for a few known NSXBs, the 
clearest example of which is the LMC transient X0535- 
668 (A0538- 66) with a peak luminosity of Lx ~ 1.2 x 
10^^ erg s"^ iBradt fc M c C lintock 198a)- 

A definitive signature of a NSXB is a detection of a 
Type I X-ray burst, a result of the accumulation of mat- 
ter onto the surface of a neutron star where it undergoes a 
thermonuclear fiash. These are not seen in black hole sys- 
tems due to the lac k of a physical surface for matter to 
accumulate onto fe.g. iNaravan fc HevJl20o3) . The duration 
of these fiashes range from seconds to minutes, with in- 
tervals on time scales of hours or days. Their X-ray spec- 
tral shape is well described by a BB spectrum which soft- 
ens as it decays as a result of the cooling of the neu- 
tron star photosphere. However, although the peak of the 
emission in strong bursts can be up to 10'^ times higher 
than t he persistent fiux. it cannot exceed the Ed dington 
limit iLewin. van Paradiis fc van den Heuvell Il995ll . Most 
importantly, NSXBs accreting at very high rates {Lace ^ 
0.25I/Bdd) are expected to be in a stable regime whe re no 
thermonuclear bursts are seen llNaravan fc Hevlll2003l) . and 
hence the absence of bursts in our source lightcurves is ex- 
pected and does not assist in distinguishing between neutron 
star and black hole primaries. 

The X-ray spectral shapes of known NSXBs and BHXBs 
may provide clues to the type of compact object present 
in these systems in MlOl. The observed spectra of bright 
NSXBs {Lx ~ 10^* erg s~^) consist of two components: BB 
emission thought to arise in the neutron star envelope with 

^ This also implies that in most cases there is no strong re- 
quirement for additional absorbing material in the vicinity of the 
sources above that through the disc of MlOl. 



kT ~ 2keV jMitsuda et al.lll984^ . and a softer DISKBB 
component from the ac cretion disc {Tin '-^l.SkcV) that 
cools with decreasing Lx l|Tanaka fc Shibazaki.l99a1 . When 
Lx < 10^^ erg s~^, the spectra become PL dominated. 
For example, the LMC transient A0538-66 displayed a BB 
spectrum with kT ~ 2.4 keV during an Einstein observa- 
tion o f the source in outburst jPonman. Skinner, fc Bedford! 
Il984h . Alternatively, if the compact object is a black hole, 
the spectral shape is characterized by a soft DISKBB com- 
ponent with lower typical temperatures (^1.2 keV) plus a 
hard powerlaw tail. Since all of our sources have Lx > 
10"^* erg s~^, we would expect any spectra of NSXBs to have 
both BB-I-DISKBB components, but this is not the case. 
The only single-component source that is well fit with a BB 
model is XMM-7, but this displays a soft spectrum with 
kT ^0.14 keV, far too cool for a neutron star envelope. To 
investigate whether the spectra of any of the two-component 
sources are consistent with the compact object being a neu- 
tron star, we have attempted to fit the data with a model 
where the cool component is represented by a DISKBB and 
the hard component by a BB. This generally gave reasonable 
fits {xt ~0.9-1.41), but resulted in DISKBB temperatures of 
Ti„ ~0. 2-0.4 keV and BB temperatures of kT ~0.1-1.5keV, 
both too low to be consistent with neutron stars accret- 
ing at the Eddington limit. The only exception to this was 
XMM-13, where the fit yielded an unrealistically high BB 
temperature of kT ~4.3keV. Hence the lack of these spec- 
tral signatures in any of the fourteen sources suggests that 
these systems are more likely to contain black holes rather 
than neutron stars. 

T.3 Single Component Spectral Fit Sources 

7.3.1 High Temperature Disc Black Body Sources 

Four of the five XMM sources that are best fit with a simple 
absorbed DISKBB model have hot inner disc temperatures 
(7in=0.96-1.33 keV), two above and two below the ULX lu- 
minosity boundary. It is now clear from RXTE studies is 
that the temperatures of these sources are clearly within 
the range observed from Galactic BHXBs in the HS (disc- 
dominated) emission state^ whic h lie in the range 0.7-2.2 keV 
l|McChntock fc R,emiUarJl2003t l. This again suggests a sim- 
ilarity between our sources and the known Galactic BHXB 
population. 

7.3.2 Supersoft Sources 

XMM-7 has a very soft blackbody spectrum with a temper- 
ature of T„,=0.14lg;°^keV (BBODY) / r,„=0.16lHl keV 
(DISKBB). Such supersoft sources (SSSs) have been de- 
tected in the Milky Way (e.g. iGrc incr 2003) as well as a 
number of other galaxi es fMlOl. IPence etal] l200ll: M81, 
ISwartz et alJl200l M 31.'Di Stefano et alj|20qa Th e Anten- 
nae. lFabbiano et all 2003b; M83. ISoria fc Wj2003l) with lu- 
minosities of lO^'^ — lO''^ erg s~^. Although SSSs were clas- 
sically d efined as sourc es with emission dominantly below 
0.5 keV jGreineJlioOoll . this has more recently been ex- 
panded to include sources with either blackbody type spec- 
tra with kT < 175 eV or powerlaw spectra with F ^ 3.5, 
with les s than 10 percent of the total luminosity above 
1.5 keV jDi Stefano fc Kond l2003bl) . This class of X-ray 
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source may include several object types, including SNRs 
and accreting neutron stars. However, it is generally be- 
lieved that most are a result of nuclear burning on the sur- 
face of white dwarfs (WDs) accreting at high rates, with 
a maximum temperature of ~150eV corresponding to the 
Eddington luminosity for nu clear burning on a IAMq WD 
JPi Stefano fc Kond l2003b^ It is therefore possible that 
XMM-7 with Lx ~ 10^^ erg s~^ is an extreme example of 
a WD SSS, but another interesting possibility is that the 
emission in extremely luminous SSSs such as this arises in 
the accretion discs of IMBHs, since such high luminosities 
coupled with such low disc temperatures imply black hole 

masses of ^ lO OMpi . 

iKing fc Poundj (|20o3) show that ultra- 
soft components in super-Eddington sources may result 
from Compton-thick outflowing material from stellar-mass 
black holes, with effective photospheric sizes of a few 
Schwarz child radius Rs, c omparable to that expected in 
IMBHs. iMukai at all 12003) use this outflow scenario to ex- 
plain the soft X-ray emission in the MlOl transient source 
P98 (not detect e d in th e JTMM observation, see section 1731 . 
iFabbiano eraP J2003d) discuss a highly variable supersoft 
ULX in the Antennae, and demonstrate that its changing 
spectral characteristics can be explained a WMq black hole 
with a mass outflow rate consistent with that expected from 
thermal timescale mass transfer in a massive XRB or for a 
bright soft X-ray transient (SX T) burst, the two scenarios 
put forward to explain ULXs bv lKine| ll2002h . 

Although we do not have sufficient data to test whether 
the spectral shape of XMM-7 varies with time, its 0.5-2 keV 
emission does show some variability on long timescales, 
which argues against the alternative interpretation for this 
source of a HNR (see section 16.71 . It has an intrinsic lu- 
minosity ~ lO'^^ erg s~^ in the XMM observation and may 
indeed be a stellar-mass XRB with an outflow. Future obser- 
vations of this source with XMM-Newton or Chandra may 
give additional clues to its nature. 



14). Although these are marginally better fit with the PL 
model (see comparison of ratios in Figure El, the DISKBB 
fits for these sources all give ^ 1 with disc temperatures 
in the range of 0.7-1.5 keV, similar to the sources discussed 
in section 17. 3. II It is therefore entirely possible that eight out 
of nine single component sources are XRBs in the HS state. 
We explore this idea further by investigating whether the 
underlying hard continuum in the multi-component sources 
can be modelled with a DISKBB in section VTM 

An extremely interesting result is that only one source, 
the nuclear source XMM-10, has an unambiguous powerlaw 
spectral shape with F = 2.3. This marks XMM-10 as an 
unique source in the sample, and could be interpreted as ev- 
idence that the nucleus of MlOl harbours a low- luminosity 
AGN (LLAGN), although the spectrum is steeper than the 
standard AGN photon index (F ~1.7-2.1,|Nandra_& Pount^ 
Il994h . XMM-10 is resolved into two sources in the Chan- 
dra data, both with powerlaw spectra (see section 10.1011 . 
with the source coinciding with the optical nucleus (within 
~35pc) possessing the steepest spectrum of F = 2.1. The 
lightcurve data for XMM-10 shows that at least one of the 
two Chandra sources is variable on both long- and short- 
timescales, but it is impossible to distinguish which of the 
two may be dominatin g the emission i n the XMM data. 

As a comparison, 'Baga noff et alJ lf2003) have recently 
reported that the Chandra spectrum of the X-ray emission 
from the Sagittarius A* massive dark object at the centre 
of our own Galaxy (morphologically similar to MlOl) can 
be fitted with a steep powerlaw (F = 2.7) spectrum, albeit 
with an extremely low luminosity of 2.4 x 10^^ erg s~^ (0.5- 
7 keV) , consistent with a quiescent super- massive black hole 
(SMBH). This is a factor ~ 10^ less luminous than XMM-10. 
Interestingly, there is no evidence for an AGN in MlOl in the 
optical regime, with the nucleus c lassified as an HII region by 
IHo. Filippenko. fc Sargeli^ (^93) based on its emission-line 
ratios. It is therefore perfectly plausible that we are simply 
detecting a bright XRB in the vicinity of the nucleus. 



7.3.3 Powerlaw Sources ? 

Another result of this analysis is that the spectra of four 
sources (XMM-6, 10, 12 fc 14) are best-fit with sim- 
ple absorbed powerlaw continua with F ~1.8-2.6. Lu- 
minous sources with powerlaw spectra (as opposed to 
DISKBB spectra) have been found in recent studies e.g. th e 
XMM-Newton minisurvey of ULXs jFoschini et al.ll2002l) : 
NGC 5204 X-1 llRoberts et aJl200ll) : M 51 (NGC 5194 #3 7 
fc 82. iTerashima fc Wilsonll2003l) : M 83 JSoria fc Wj2003l) : 
NGC 4485/90 iRoberts et al.ll20o3) . While powerlaw spec- 
tra with F=1.5-2 are typical of the low/hard (LH) state of 
Galactic BHXBs llTanaka fc Shibazaki|[l99^ . steeper spec- 
tral shapes could represent strong Comptonized thermal 
emission from a black hole accretion disc in a very high/steep 
powerlaw state (VHS/SPL). Ba sed on RXTE observations, 
iMcChntock fc RemiUardI (l2003h define the SPL state by the 
presence of a powerlaw component with F > 2.4. Out of 
these four XMM sources, only XMM-6 & 10 fall within this 
category (within the 90 percent confidence limits), while 
XMM-12 fc 14 have harder powerlaw slopes of F ~ 1.8 - 1.9 
consistent with the LH state. 

However, in this dataset we cannot reject the DISKBB 
model fits for three of four of these sources (XMM-6, 12 fc 



7.4 Multi-component Sources 

Five of our luminous sources are best fit with multi- 
component spectral models. Similarly, two-component mod- 
els have been required to adequately fit some other re- 
ce iit XMM and Cha ndra data e.g. NG C 1313 X-1 and X- 
2 jMiUer et al.ll2003l) : NGC 6 946 X-11 llRoberts fc Colberli 
l2003h : two M51 ULXs fTerashi ma fc WilsonI l20ol) : 
NGC 5408 X-1 jKaaret et al. 2003). Even though our most 
luminous source (XMM-1) is well fit with a simple DISKBB 
model, the other sources provide evidence that many sources 
could have more complicated spectral shapes than are shown 
by the simple fits that can be done on photon-limited data. 

To begin with, we fitted these sources with the canon- 
ical model of a powerlaw continuum plus a soft compo- 
nent to model the soft excess. Three are best fit with 
PL+DISKBB/BB models (XMM-2, 11 & 13), although 
there is little difference statistically between the DISKBB 
and BB fits. The disc temperatures are cool, and in two cases 
(XMM- 11 fc 13) are similar to the inner disc tempera tures 
of SSSs discussed in section ITX^ iMiller et al.l (l2003l) find 
similar disc temperatures in their study of two ULX sources 
in NGC 1313, and interpret this as evidence of IMBHs based 
on the argument that the black hole mass scales inversely 
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with the accretion disc temperature (Ti„ cx M~^^'^). The 
ambiguity of which type of blackbody model is the best fit 
plus the fact that our source luminosities do not imply IMBH 
masses casts some doubt on this interpretation in our data. 
Interestingly the powerlaw slopes are all < 2 (XMM-13 has 
a particularly hard slope with r=0.36), harder than seen 
in the HS state of Gal actic BHXB candidates (r=2. 1-4.8, 
iMcClintock fc Remillard 2003) , and all have a substantial 
contribution to the total source fluxes. 

The other two sources (XMM-3 & 5), both with ULX 
luminosities, are best fit with PL+MEKAL models. As dis- 
cussed in sections 16.31 & 16.51 both are located far from 
the central region of the galaxy where most of the dif- 
fuse thermal gas is found, and so it is unlikely that what 
we detect here is contamination from the diffuse emis- 
sion in the disc of MlOl. We are therefore likely to be 
detecting thermal plasma in the vicinity of the hard X- 
ray source. Although this could suggest the presence of 
SNRs, none was d etecte d at the location of XMM-3 by 
iMatonick fc FesenI lll997ll . though the position of XMM-5 
is just outside the field covered in that study. In the case 
of XMM-3, the MEKAL temperature is hot (fcT=0.98keV) 
and it may be the case that we are detecting hot photoion- 
ized plasma surrounding a high-mass companion star in a 
binary system. Galactic high-mass XRBs (HMXBs) can ex- 
ihibit strong recombination lines if accreting from hi ghly 
ionized stellar wind s (e.g. Cyg X-3. ISo ria fc Wull2003l and 
references therein). iTerashima fc w'iisonril200^r point out 
that the luminosity of this thermal component is likely 
to constitute < 10 percent of the total flux of the sys- 
tem, which is consistent with the 5 percent MEKAL con- 
tribution we measure here. In the case of XMM-5, the 
plasma component makes a much more substantial contri- 
bution to the total source flux (20 percent), but its temper- 
ature is much lower (fcT=0.17keV), and we could be see- 
ing thermal emission related to star formation reminiscent 
of the low-temp erature component of starburst superwinds 
(e.g. NGC 253. 'Pietsch et al.l l200ll IStrickland et alJl2000l: 
NGC 3256. iLira et al. 200211 . although lack of evidence for 
star formation activity at other wavelengths is problematic. 
The powerlaw components in these sources are soft (F > 2), 
similar to those of XMM-6 fc 10 discussed in the previous 
section. XMM-13, located in the giant HII region NGC 5461, 
is another interesting case, as mentioned in section lli.131 Al- 
though not fully constrained with these data, the best fit is 
a two-temperature MEKAL plus hard PL model (F=0.3, 
fcT=0.3/1.0keV), where the MEKAL components consti- 
tute ~43 percent of the total fiux of the source, again rem- 
iniscent of starburst superwind plasma temperatures which 
could perhaps arise in a very energetic HII region. 

Since we are speculating that all of the single- 
component spectral fit sources in MlOl (except the nuclear 
source XMM-IO) are accreting XRBs in a high state (see 
section 17. 3. 31 . we have also attempted to model the un- 
derlying hard continuum in the multi-component sources 
with a DISKBB with temperatures similar to those found 
in stellar-mass black holes (~l-2keV). In three cases 
(XMM-3, 5 fc 11), this was successful, with statistical im- 
provements at the > 99 percent level compared with the 
single-component DISKBB fits. XMM-3 could be modelled 
by both DISKBB-I-PL (Ti„=0.9keV, F=3.2, xl ~ 1-1) 
and DISKBB-FBBODY (Ti„=0.9keV, fcr=0.2keV, xl ~ 



1.0) models; XMM-5 with a DISKBB-f MEKAL model 
(Ti„ = 1.2, fcT=0.2keV, xl ~ 1-4) and XMM-11 with a 
DISKBB-HBBODY model (r,„=1.5keV, fcT=0.2keV, xl ~ 
1.3). The fits for the remaining two sources (XMM-2 fc 
13) yielded unrealistically high disc temperatures (> 3keV). 
This interpretation of the data has important implications 
for the IMBH scenario. The hard continuum in such sources 
is usually interpreted as powerlaw emission similar to that 
seen in the LH state of Galactic XRBs; if the soft excess 
is modelled by a cool DISKBB, the best-fit temperatures 
imply black hole masses in the IMBH range. However, if 
the hard component can be equally successfully modelled 
by a DISKBB with a "normal" stellar-mass black hole tem- 
perature, then this component could represent the thermal 
signature of the accretion disc (with a temperature typical 
of a stellar-mass black hole) and the soft component must 
have a separate origin. Again this implies that an IMBH 
is not required. Further detailed spectral studies of similar 
sources are required to explore this possibility. 

T.5 Transient Sources 

It is believed that most low-mass XRBs (LMXBs) will ex- 
hibit transient behaviour during their lifetimes (King 2002j), 
in the form of SXT outbursts caused by hydrogen-ionization 
instabilities in the accretion disc around the black hole or 
neutron star. The most obvious transient source we have 
detected in this study is XMM-2, which has been in a qui- 
escent state in all previous observations and has brightened 
by a factor of ~30 in the XMM observation. It also is worth 
commenting on the lack of a detection in the XMM observa- 
tion, and therefore the transient behaviour, of the most lumi- 
nous ULX (P98) in the Chandra observation, a highly vari- 
able source which seemed to reflect anticorrelated ch anges 
in the source temperature and size jMukai et al.ll2003h . As- 
suming the model parameters of lPence^^Lril200 j) . we cal- 
culate an upper flux limit of ~ 8 x 10~^^ erg cm~^ s~^ 
in the XMM data, which implies a drop in luminosity of 
at least a factor of ~20 since the Chandra observation. 
Such transient sources are likely to be LMXBs, since known 
black hole HMXB systems are persistent X-ray sources 
(McClintock fc Remillard 2003)®, although the caveat here 
is that the known black hole HMXB systems are wind ac- 
cretors, i.e. they are not going through a phase of thermal- 
timescale mass t ransfer that may be ne cessary to reach ULX 
luminosities (see iKaloeera et al.ll2003l for more detailed dis- 
cussions). Future monitoring observations of galaxies such 
a MlOl, designed to detect transience over timescales of 
months, might therefore provide a means of distinguish- 
ing between LMXB and HMXBs possessing black hole pri- 
maries, which cannot be distinguished spectroscopically in 
the X-ray regime. 



8 CONCLUSIONS 

In this paper, we have studied the X-ray spectral and tem- 
poral properties of the most luminous sources in MlOl at 

^ Note that though SXT outbursts may occur in HMXB systems 
with Be- star companions, these systems possess a neutron star 
primary iTauris fc van den Heuve^l2003^ . 
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the time of the XMM-Newton observation, and our results 
suggest that they are all accreting systems associated with 
MlOl. For comparison purposes, we have also analysed two 
archival Chandra observations of MlOl to investigate spec- 
tral and temporal changes in these sources. The positions of 
thirteen of the fourteen XMM sources are covered by either 
or both of the Chandra observations. Of these, eleven are 
unresolved by Chandra, while the remaining two (XMM-2 
and the nuclear source XMM- 10) are both resolved into two 
discrete sources. Our results can be summarized as follows: 

• Fourteen sources in the XMM observation have suffi- 
cient counts for spectral fitting, with intrinsic luminosities in 
the range of 3.4x 10^**-3.4x 10^^ erg s~\ If these are BHXB 
systems accreting at the Eddington limit, these luminosities 
imply black hole masses of ~ 2 — 23Mq , consistent with the 
measured range of Galactic stellar-mass black holes. 

• The locations of the sources suggest that they are all 
associated with MlOl (HII regions/spiral arm/nuclear), and 
suggest a link between the high-luminosity XRB population 
and young stellar populations. 

• Five sources are best described by simple DISKBB 
models. Four have hot inner disc temperatures 
(Tiri=0.96-1.33 keV), consistent with the disc temper- 
atures found in Galactic BHXBs (Ti„=0.7-2.2 keV, 
iMcChntock fc R,emillardll2003h in the HS state. 

• Four sources have simple absorbed powerlaw spectra 
with r ~1.8-2.6, two of which are consistent with the SPL 
state. However, three of these sources are also well modelled 
by DISKBB models with disc temperatures of 0.7-1.5 keV, 
making it entirely posstble that eight out of nine sources 
single-component sources are HS XRBs. Only the nuclear 
source (XMM- 10) has an unambiguous powerlaw spectrum, 
which may be evidence of either an LLAGN or alternatively 
an XRB (given no indication of AGN activity at other wave- 
lengths) . 

• Five sources require multi-component spectral fits, in 
the form of an underlying hard PL continuum plus a soft 
excess. In three sources, the soft component can be mod- 
elled by either DISKBB/BBODY models, two of which have 
cool disc temperatures which could be interpreted either as 
evidence for the presence of IMBHs or as the signature of 
a Compton-thick outflow from a stellar-mass black hole. In 
two (possibly three) sources, the soft excess is well modelled 
with one or more MEKAL thermal plasma components. The 
hot (fcr=0.98keV) component in XMM-3 could arise in a 
hot photoionised plasma surrounding a high-mass compan- 
ion star, and the cooler (A;T=0.17keV) component detected 
in XMM-5 could be thermal emission related to star for- 
mation activity with a temperature similar to those of star- 
burst superwinds. If we accept the PL-|-(2xMEKAL) inter- 
pretation of the data for XMM-13, we could be detecting a 
multi-temperature plasma from the superbubble encompass- 
ing an XRB. In three cases (XMM-3, 5 & 11), the underlying 
hard component can also be successfully modelled by a hot 
DISKBB with a "normal" stellar-mass black hole accretion 
disc temperature (Ti„=0.9-1.5), arguing against the IMBH 
scenario. 

• Eleven sources show statistical evidence of short-term 
variability during the Jf'MM observation. The long-term data 
demonstrate that most are also variable over a baseline of 
11-24 years, with the majority varying between a factor of 



^2-4 in observed luminosity. This variability is strong evi- 
dence in favour or these sources being accreting XRB sys- 
tems. 

• In the absence of the spectral signatures expected from 
neutron stars in outburst, we conclude that the accreting 
objects in these binary systems are more likely to be stellar- 
mass black holes. 

• We find no evidence to support the HNR scenario in 
those sources which have previously been suggested to be. 
Evidence of long- and short-term variability coupled with 
the DISKBB spectral shape of XMM-1 effectively rules out 
the HNR scenario for this source. Although we find no con- 
clusive evidence to rule out XMM-7 (NGC 5471B) as a HNR, 
it has a supersoft XMM spectrum, which may result from 
nuclear burning on the surface of a WD accreting at a high 
rate, a Compton-thick outflow of material from a stellar- 
mass black hole or the accretion disc of an IMBH. 

• Even though the sources show a variety of spectral 
shapes and hardness, there is no apparent spectral distinction 
between those above and below the ULX luminosity thresh- 
old (10"^^ erg s~^), implying that we are seeing the extreme 
high-luminosity end of a "normal" XRB population. 

• Out of the nine sources with spectral fits in both the 
XMM and Chandra observations, no major spectral changes 
are apparent. While the majority do show some degree 
of spectral softening with increasing luminosity as seen in 
Galactic BHXBs, the transient source XMM-2 shows little 
change of hardness despite a significant change in luminos- 
ity. 

• Two transient sources have been detected in these ob- 
servations; XMM-2 in the XMM observa t ion an d the Chan- 
dra source P98 studied bv iMukai et alJ ll2003l) which then 
fell below the detection threshold of the XMM data. We 
speculate that since LMXBs are known to exhibit transient 
behaviour, future multiple observations of galaxies such as 
MlOl may provide a way to distinguish between LMXBs 
and HMXBs. 

We have now reached the era of X-ray astronomy in which 
detailed studies of the most luminous individual XRBs be- 
yond the Local Group (SMC, LMC, M31) are possible. This 
paper demonstrates that these sources, the extreme exam- 
ples of accreting stellar remnants, display a wide range of 
X-ray characteristics. With the superb capabilities of XMM- 
Newton and Chandra, we now have an exciting opportunity 
to begin to expand our knowledge of these and other objects 
apparently accreting at or above the Eddington Limit. 
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